High order quantum decoherence via multi-particle amplitude for boson system 
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In this paper we depict the high order quantum coherence of a boson system by using the 
multi-particle wave amplitude, whose norm square is just the high order correlation function. This 
multi-time amplitude can be shown to be a superposition of several "multi-particle paths" . When the 
environment or a apparatus entangles with them to form a generalized "which- way" measurement 
for many particle system, the quantum decoherence happens in the high order case dynamically. An 
explicit illustration is also given with an intracavity system of two modes interacting with a moving 
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I. INTRODUCTION 



| A most profound concept in quantum mechanics is coherent superposition of quantum states. This is obviously 
^vq , reflected by the special interference of two or more "paths" in terms of single particle wave function. However, this 
t-H ' coherence phenomenon does not sound very marvelous for the same circumstances can also occur in classical case, such 
as an optical interference. To manifest the intrinsically quantum features of coherence beyond the classical analogue, 
Glauber's higher order quantum correlation function (HOQCF) was introduced [1] to accounts for such higher order 
quantum coherence effect in the Hanbury-Brown-Twiss experiment [2] . 

The quantum coherence both in the first order version and its higher order one shows the embodying of the wave 
*G ' nature in the world of the microscopic particles, but it is not robust due to the wave function collapse or quantum 
decoherence caused by a quantum measurement or by the coupling environment. Roughly, in the single particle picture, 
' this decoherence phenomenon losing coherence can be understood according to the quantum entanglement of the 
considered system with the environment or the measuring apparatus. Obviously, this entanglement implies a "which- 
path" detection[3]. Precisely speaking, in a initial coherent superposition \tp s ) = J2c n \n) , each system state \n) 
q-( corresponding a "path" is correlated with an environment state |e„) to form an entangling state \4>t) = X) c n\ n ) & l e «)- 
• • ■ The existence of the different states |e„) distinguish among the "paths" of different \n). The interferences terms in 
I(x) = Tr((x\?pT) (iPt\%)) can disappear when the environment state \e n ) are very distinct, i.e., each path is explicitly 
labelled (e m |e„) = 5 mn . 

For the first order decoherence, the above well-known explanation in terms of "wich-path" detection mechanism is 
more simple, but very profound. However, it is not obvious yet whether this mechanism can account for decoherence 
in the higher effects: This is because we do not exactly know what is the "paths" and the corresponding "which-path" 
detection in the high order version. More recently, we have touched the quantum decoherence problem in higher 
order case [4]. The concrete calculation motivated us to further consider the "which-path" picture for the higher 
order quantum decoherence. In this paper, for a boson system we introduce the concept of the multi-particle wave 
amplitude, whose norm square is just the high order correlation function. As an effective wave function, this multi-time 
amplitude can be shown to be a supposition of several components. When the environment or a apparatus entangles 
with them, the quantum decoherence occurs in the high order case dynamically. This decoherence process losing the 
higher order coherence can be explained as a generalized "which-path" measurement for the defined multi-particle 
paths. 



II. MULTI-PARTICLE PATHS FOR PROBABILITY AMPLITUDES 



Firstly, we start with considering the meaning of the "path" in the high order quantum correlation. The typical 
example of the higher order quantum coherence [5] is that the single-component state |lfc, of the two independent 
photons with momenta k and k! shows its quantum coherence in its second order quantum correlation function 



G^ = G^ 2 \r\,r 2 ,ti,t 2 ). It can just be written as of form of the norm square G^ = 2 of the equivalent "two-time 
wave function" 

^ = i>(r u r 2 ,h,t 2 ) = (00|£^(r 2 ,t2)E + (ri,ti)|l fc ,lfc>> (1) 

It was also called the biphoton wave packet for the photon field E + (r, t) [6] . Especially, we remark that ip is a coherent 
superposition of several "probability amplitudes" . This result can be promoted to the universal case with quantum 
systems of identical particles. 

To see the main physical ideas implied by the above well-known result, without loss of the generality, we define a 
"measuring" operator of two modes V and H [4] 

4> = c v b v e-^ vt + c H b H e- iuJHt = c v {t)b v + c H {t)b H (2) 

where bu and by are the annihilation operators of the boson system. The generalized second order correlation 
function[5] 

= | <0, O|0(t 2 )^(ti) 1 1 v, Iff > | 2 = |*(ti,*2)| 2 (3) 

The two time wave function &(ti,t 2 ) can be understood in terms of the two "paths" picture from the initial state 
to the finial one |0, 0) : 





\i v ,0h) c ^ 2) |o,o) 




c H (t 2 ) y 
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They are just associated with the two amplitudes forming a coherent superposition 

*(ti, t 2 ) = c v c H e- luJvt2 - tuJHtl + CHCve-^H^-^v^ (4) 

Correspondingly, the second order correlation function 

G (2) = 2\c v c H \ 2 [l + cos(K - u H ][t2 - ti])] (5) 

The above observation for the second order quantum coherence can also be discovered in the higher order case. 
Our arguments in this paper will be based on two novel observations: a. The equivalent field operator <l = ^ c n b n is 
specified for a quantum measurement to a superposition state \<f>) = ^2 c n \n). b. For a certain initial single component 
state I so) of N particles system, the n'th order quantum correlation function 

G^(r 1 ,r 2 ,---,r n ,t 1 ,t 2 ,---,t n ) = \^\ 2 (6) 

can be written as the norm square of an effective wave function ip( n ') , which is just a superposition of many amplitudes. 

For the seek of simplicity, we only consider the third order situation when the initial state is in the state 2y). 
Indeed, the generalized third order correlation function 

- (2 v lH\$Ht 1 )$\t 2 )^\t 3 )^(t 3 )$(t 2 )$(t 1 )\2 v ,lH) 

= |(0,0|^3)^2)«A(il)|2y,lH)| 2 = |f(ti,t2,*3)| 2 (7) 

is indeed a norm square of the two time wave function : 

*(ti,t2,*3) = V2c v 2 c H e- luJv{t3+t2) - WHtl + 

V2c H c v 2 e- luJHt2 - luJv( - t3+tl) + y/2c H c v 2 e~ lulHt3 ~ lu;v{t2+tl) (8) 

Each term in the above equivalent wavefunction is contributed by the corresponding one of the four "paths" from 
\2vAh) into |0,0): 
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From the above equivalent 3-time-wave function, the third order correlation function is explicitly written down 

= A\c v 2 c H | 2 (^ + cos[(^y - iu H )(t 2 - ti)] + (9) 

C0s[(w y - U) H )(t 3 - il)] + COs[(^y - U) H )(t 2 - t 3 )D 

and shows the quantum interference in the time-domain. 

From the above calculations for the second and third order quantum decoherence, the observation can be made 
that, for a specially-given initial states, a high order correlation function may be explicitly written down as the norm 
square of the equivalent multi-time-wave function, which is a coherent superposition of several complex components 
in associated with the generalized many-particle paths. It is pointed out that this kind of many-particle path is not 
a simple-product of single-particle paths and it can be determined by the concrete measurement. 

III. INTRACAVITY MODEL 

With the above introduced concept of "many-particle path" , we can discuss the higher order decoherence problem 
by considering the many-particle "which -path" measurement. Let us use the following model to sketch this central 
idea. Our model is consisted of non-dissipative bosons in two modes. The problem is studied via calculating the 
second order quantum correlation functions in the Heisenberg picture. We take h = 1 in this paper. 

The model Hamiltonian H = Ho + V is defined by 

Ho = uybybv, 

V = y^Uja\aj + y~l[d v (u)j)b\rb v + dH(u)j)b\jb H ]{a] + dj), (10) 

3 j 

where Ho is the free Hamiltonian of the system, V the free Hamiltonian of the reservoir (or a detector) plus the 
interaction between the system and the reservoir; and b v (t>v), V H (bu ) the creation (annihilation) operators for two 
modes with frequencies ujy an d &h = . The operators aj(aj) are creation (annihilation) operators of the reservoir 
modes of frequencies u>j . The frequency-dependent constant dui^j) (dv(wj)) measures the coupling constant between 
H(V )mode and the j mode of the reservoir. The most important feature of the model is that [JTo,^] = 0, i.e. the 
system does not dissipate energy to the reservoir, but it can leave imprinting on the reservoir since, for different 
number states |ny,n#), there are different interactions ji n v dv 'j) + nndn (o>j)](aj + dj) acting on the oscillator 
reservoir with different driving forces <~ nydvi^j) + nHdni^j)- 

When there is only one mode in external system we can physically realized this model as an intracavity model[4,8]: 
Two mode cavity field interact with a moving wall of the cavity, which is attached to a spring and can be regarded 
as a harmonic oscillator with a small mass. The coupling of fields to the cavity wall (a moving mirror) is just by the 
radiation pressure forces proportioned to the photon numbers by V H i>H and by'b v . 

The second order coherence is directly determined by the second order correlation. 

G[t, t', (3(0)] = Tr(p(0)B\t)B\t')B(t')B(t)), (11) 

which is defined as a functional of the density operator ps(0) of the whole system at a given time 0. Here, the bosonic 
field operator 

B(t) = exp(iHt)[ab H + c 2 b v ] exp(-iHt) = (12) 

exp(iVt)[cibH + c 2 bv exp(— iuyt)] cxp(-iVt) (13) 

describes a specific quantum measurement [4] with respect the polarized photon states |+) = Ci|i?) +C2IV") and |— ) = 
C 2\H) — C\\V) where ci and c 2 satisfy the normalization relation |ci| 2 + |c 2 | 2 = 1. Without loss of the generality, we 
take c\= c 2 = 1/V2 standing for a given measurement as follows. 
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To examine whether the macroscopic feature of the reservoir causes the second order decoherence or not , we 
consider the whole system in an initial state 

|V(0)) = |1h,M®|{0,}), (14) 

where \{0j}} is the vacuum state of the reservoir. Here, we have denoted the general Fock states of the many mode 
field by \{nj}) = \ni,n 2 , ...). 

For the total system, in stead of defining the equivalent "two-time wave function" in the above section , we define 
an effective two-time state vector 

\^ B (t,t'))=B(t')B(t)\m)- (15) 
to re-write the second order correlation function as 

G[t,t',p(0)] = ^ B (t,t'MB(t,t')) (16) 
It is interested that the effective state vector can be evaluated as the superposition 

|V> S (M')> - ie^( -°) t '[exp(-^ y i')e- l,>(1 ' 0)t 'e^ (1 ' 0)t + (17) 

exp(- i ^ y t) e ^(°' o ) t ' e - l,> ( o ' 1 ) t ' e ^( o ' 1 ) t ]e-^( 1 ' 1 ) t |{0 J }) ® \Q H , Oy) (18) 

of two components for the two paths from the initial two particle state \1h, ly) to the two particle vacuum \0h, Oy) 
. It should be noticed that the effective actions of the reservoir 

V(m,n) = Vj(m,n) — V^Wjajaj + V\dy(wj)ra + rfjj(a;j)n)(at + hj) (19) 

3 3 j 

can label the different paths and thus lead to the higher order quantum decoherence. The above result 
clearly demonstrates that, in presence of the reservoir, the different probability amplitudes (~ cxp(— iujyt') and 
exp(-iwyi)) from \1 H , ly) to |0 H ,0y) entangle with the different states ( \ e iV(o,o)t' e -iV(i,o)t' e iV(i,o)t e -iV(i,i)t^ Qj ^ 
and i e ^(o,o)t' e -iy(o,i)t' e ^(o,i)t e -iV(i,i)t|{o j .}) ) G f the reservoir. This is just physical source of the higher order 
quantum decoherence. In the following section an explicit calculation of the second order correlation function will be 
given to show this crucial observation. 



IV. DYNAMIC DECOHERENCE IN HIGHER ORDER CASE 

In our calculation, the second order correlation function 

G[t, t>, p(0)] = \ + i e --(t-0 Y[ Fj + \e~^ (*"*') J] F*, (20) 

3 ' 3 

is firstly expressed as a factorization form [9] where each factor 

Fj = (0,|^(t 6 )|0,) = (o^a^V^ ' 1 )^ '^ (21) 

is a two-time transition amplitude of the j'th mode of the reservoir. Obviously, the term f\- Fj determines the extent 
of coherence and decoherence in the second order case, which is called "the decoherence factor" and plays the role 
just as the same as that in the first order decoherence [9]. 

In the following, to given the factor Fj explicitly , we adopt the Wei-Norman method [10,11] to calculate the 
equivalent time evolution defined by u 6 j{t§). It can be image as an evolution governed by a discrete time-dependent 

Hamiltonian H{t) dominated by Vj(l, 1), —Vj(l, 0), Vj(l, 0), — V}(0, 1), Vj(0, 1) and — Vj(l, l)in six time-intervals [to = 
0, h=t}, [t!,t 2 = 2t], [t 2 , h = 2t + t'), [t 3 ,t 4 = 2t + 2t'], [t 4 , h = 3t + 2t'], [t 5 ,t 6 = At + 2t'\ respectively. In the fc-th 
step of calculation, we take the final state of (k — l)-th step as its initial state. Therefore,we obtain u®(te) as the sixth 
step evolution 
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Here, g^.j(te)(k — 1,2,3,4) are the coefficients that can be explicitly obtained, but for the calculation of the j-th 
component Fj — exp[glj(t 6 )] of the decoherence factor, we only need to know g\y The detailed discussion in the 
appendix gives 



94j{t6) = -o sin [ j + l -2 

3 3 

[uj(t' -t) + 2(1 - cos(w j [i / - t]) sinwji] + (1 -2cosUjt) sin^i' - t})]. 



(23) 



Since the real part of g| ■ is no more than zero, the norm of Fj is no more than one. From the same argument as in 
the first order decoherence [9] , the universal factorization structure of the decoherence factor implies the second order 
decoherence in the macroscopic limit. 

In order to demonstrate the above arguments quantitatively, we give the numerical results for the second order 
decoherence for different numbers N of the quantum oscillators. As N increase, these results are illustrated in FIG.l. 
In the numerical calculation, the coupling constants {dv(j)} take random values in the domain [0.8, 1.0], the coupling 
constants {dn(j)} in [0.2,0.4], and the frequencies {uij} in [0.5, 1.5]. The other parameters is given in the caption of 
the figure. 
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FIG. 1. The horizontal axe denotes time period t' — t, the vertical axe denotes the second order correlation function 
G[t, t', p(0)], parameters cu v = 1.0, (a)N = 5, t = 0, (b)N = 5, t = 5, (c)N = 10,t = 0, (d)JV = 10, t = 5. 

From the above illustration, we see clearly that the second order correlation is not only a function of the time period 
t' — t, but also a function of time t. The amplitude of the second order correlation function is mainly determined 
by the time period, but the phase is determined by the two parameters in the same time. More important, with the 
increasing of the number of the quantum oscillators, the second order coherence vanishes faster and faster, and the 
quantum revival amplitude becomes smaller. With extrapolation, it can be expected that, when the number of the 
quantum oscillators limit to infinity, i.e., in the macroscopic limit of the reservoir, the second order coherence will 
decoherence in a short time and no quantum revival can be observed. 
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APPENDIX: 



In this appendix, the Wei-Norman method [10,11] is adopted to calculate the second order decoherence factor Fj. 
The calculation is completed in six steps. During the time period [tk-i,tk], the evolution for Wj(t, t') is dominated 
by the Hamiltonian 

h)=c§a)a i +$a)+T k j a i , {k = 1, 2, • • • , 6}. (Al) 
The coefficients a k , [3 k , 7^ and the times intervals T^=t\~— tk-i take different values in different steps: 



a) 


= Wj,0j = 


7) = 


d v {ui 3 ) + d H {uj),Ti = t, 




= ~Uj,0j 


= 7* 


= d v (uj),T 2 = t, 


«? 


= <*>,, 05 = 


7? = 


d V (Wj),T S =1?, 


a) 




= 7* 


= -d H {^j),Ti = t\ 






7f = 


d H (oJj),T 5 = t, 






7? = 


-d v (uj) - d H (vj),T 6 = t. 



(A2) 

Due to the fact that ajdj, a], dj, 1 form a closed algebra - the Heisenberg-Wely algebra, the unitary time evolution 
operator at every step takes the following form 

u) (T) = e& (T)a ' e& (T)a * ^ e4 ^ e 4 ^ . (A3) 

for T e [tfe-i,tfe] in a special sequence. The benefit of the above form is that only the coefficient g\j{T) is needed to 
be known in the calculation of the average value of the vacuum state. 

According to the Schrodinger equation i-^fU k — h k Uj, the coefficients satisfy the equations: 







d k 
dT 92 i 
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d k 
dT 9 ^ 


9ij 


d k 

dT 92j 
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"4 


d k 
dT 933 


= 


dp, k - 

jrp94j - 9ij e 


"4 


d k 
~dT 9 ^ 


= 



(A4) 



Using the results 



±gk j = - ia k g k j _ k A. g k j = _ h k gkj , • . , . 

obtained by simplifying the system of equations , we get the solution 



^•m = (^(i fe -i) + §)e-^ (T - tfc - 1 -5 



„k (rp\ _ k 1. \,rf(„k(, \ , s/ -iot k AT-tk-l) ->\ 1 „- ^7j fc ( r t k -l) 

9i 3 \T) - g 4j {tk-i) + -T{9ij{tk-i) + -r)(e > y >-l)+i -r (A6) 



Notice that we have used the step-initial conditions 

g k lj (tk-i) = 9 k lj - 1 (t k -i), (A7) 
9 k 4 j (t k -i)=9%i 1 (tk-i) (A8) 

and the initial conditions <7ij(io) = 9%{to) = 0- Then we obtain a set of iteration equations 



G 



9ij{tk) = {9i 3 [tk-i) + -feje - -r, 

3 3 

= tfoHtu - 1) + 4kir 1 (**-i + 4)(«-^ Tfc - 1) + ^4^- (A9) 

Iterating six times with different initial conditions and coefficients, the final result of g\{to) is obtained as the 
equation (23). 
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